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A UNIFIED TREATMENT OF LUNAR THEORY 


AND ARTIFICIAL SATELLITE THEORY 
By William A. Mer smart 
Ames Research Center 


SUMMARY 


Lunar theory and artificial satellite theory are treated by a unified 
method that is a generalization of the von Zeipel procedure. The technique of 
separation of variables is used to generate a single canonical transformation 
that eliminates the time and all the angle variables from the Hamiltonian. 

The validity of the method is established and the calculations are carried far 
enough to illustrate the techniques involved. 


INTRODUCTION 


The theory of artificial satellite motion has been developed in recent 
years by methods that had not previously been considered applicable to lunar 
theory. The problems appear quite different, mathematically, because of the 
explicit appearance of the time in the lunar problem. 

In the present paper the familiar von Zeipel procedure is generalized in 
two respects. Explicit dependence on the time is permitted, in a certain 
restricted form, and then the time and all the angle variables are removed 
from the Hamiltonian by means of a single canonical transformation. This 
generalization is capable of solving the complete lunar problem far more effi- 
ciently than the older methods of von Zeipel or Delaunay. It is then shown 
that this method can be applied directly, in a simplified, degenerate form, to 
the artificial satellite problem. 

The familiar Delaunay variables are used, with one modification. The 
longitude of the node is measured from a moving, rather than from an inertial 
reference direction. For lunar theory the reference direction is that of the 
sun, for artificial satellite theory, it is the meridian of Greenwich. This 
device removes the explicit time dependence from the satellite problem, and 
drastically simplifies the lunar problem. 

The generalized procedure is presented in the next section, and its 
application to lunar and satellite theory follow in that order. 



U CO 


SYMBOLS 


a 

semimajor axis 

G 

/pa(l - e 2 ) 

e 

eccentricity 

H 

G cos I 

£ 

true anomaly 

I 

inclination angle 

g 

argument of perigee 

Jn 

dimensionless coefficients in 



zonal harmonics 

h 

longitude of node from moving 


/pa" 


reference line 

L 

k 

gaussian gravitational 

M 

mass of moon 


constant 

P n 

Legendre polynomial 

l 

mean anomaly 

p m 
n 

associated Legendre function 

m 

dimensionless ratio of mean 



motions 

R 

equatorial radius of the 




earth 

n 

mean motion of moon or 




satellite 

S 

determining function, mass of 




sun 

P 

angle variable 

T n,m 

tesseral harmonic 

q 

action variable 

IV 

determining function 

r 

geocentric distance 

X,Y,Z 

decomposition types 

t 

time 

6,6 

solar factors 

u 

eccentric anomaly 

Y 

1 - cos I 

X 

action variable 

r 

longitude of sun’s perigee 

y 

angle variable 

£ 

eccentricity of sun’s orbit 

B 

earth-moon barycenter 

9 

geocentric latitude 

c n,m4 

s 1 

dimensionless coefficients in 

X 

mean anomaly of sun 

°n,m f 

tesseral harmonics 

D n 

zonal harmonic 

A 

geographic longitude 

E 

new Hamiltonian, mass of 

P 

gravitational parameter 


earth 

old Hamiltonian 

V 

mean motion of sun, spin 
velocity of earth 

F 
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p 

4 


distance from sun to earth-moon ip 
barycenter 

Q 

true anomaly of sun 


elongation of moon 

longitude of node from inertial 
reference line 


THE GENERALIZED VON ZEIPEL TRANSFORMATION 


Recent applications of the von Zeipel transformation have been restricted 
essentially to Hamiltonians in which the time does not occur explicitly 
(refs. 1 and 2). When it does occur, it is formally removed by the artifice 
of adjoining an additional pair of canonical variables (ref. 3, pp. 530-531, 
and ref. 4) . A generalized method is developed here that is applicable to 
certain problems involving a time -dependent Hamiltonian. In the case of lunar 
theory this permits many of the solar terms to be expressed in closed form, 
thereby avoiding certain Fourier series expansions that frequently complicate 
the problem unnecessarily. 

Consider the canonical system 

dXj _ 3p_ d Xi _ 3F_ 

dt 3y^ * dt 3x^ 

with M action variables” x = (xi, X 2 , X 3 ) and ’’angle variables” 
y = (y i# y y . Let the Hamiltonian be of the form 


F = F(x, y, A) 


where A (an angle) is a linear function of the time, with v = dA/dt a 
constant. This is the only form of time-dependence that will be considered. 

Now let q = (q x , q 2 , q 3 ) and p = (p 1 , p 2 , p 3 ) be new canonical variables 
defined by the implicit transformation equations 


_ 3W(q, y. A) 

p i a qi 


x - 3 W U, y. hi 

1 3 Yi 


with the "determining function," W, given by 



(Einstein's summation convention will be used throughout this paper: a 

repeated subscript is to be summed over its range.) Thus, in vector notation 


x = q + Aq 
y = p + Ap 

with the increments having components given implicitly by 

(4q) . = »c<w, » 
x) 

Expanding the right member of the last equation in a Taylor series near y = p 
gives 


Wi- as V ~ 


j ‘ 2 3q. 8 P j 8p k (Ap) j (Ap) k • • • 


This can be solved iteratively to any desired degree of accuracy, if S 
and its derivatives are assumed small: 


First order: (Ap)^ = - 


3S 

3q. 


Second order: 


Third order: 


CAp). 
(Ap) i 


3S 

3 2 S 

/ 3S 

9c k ‘ 

3qi 3p j 

\ 

3S 

3 2 S 

/ 3S 

3qi " 

3 qi 3 P j 

\ 3 <j 

1 

3 3 S 

( 3S 

2 3 qi 

3Pj 3 Pk 

V ^ 


Similarly, 


8S 


3 Z S 8S 


3 Z S 


d ^S 8S 

t ' Aq '* i ~ 3Pi 3 Pi 3 Pj 3c lj + 3 Pi 9 Pj 9c lj 9 Pk 9c (k 


3 3 S 


3S 3S 


2 3p i 3pj 3pk 3qj 3q k ' ‘ ' 
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Thus the equations 


x = q + Aq 
y = p + Ap 

with Aq, Ap given above, constitute an explicit solution of the implicit 
transformation equations. 

Now, q and p are canonical variables: 

dqi _ 3E dPi _ ^ 3E 

dt 9p^ ’ dt - 3q^ 

with the new Hamiltonian, E, given implicitly by 


E = F (x, y, X) 


3W(q, y, X) _ 3S(q, y, X) 

"at 3X 


(ref. 5, ch. 6) 

If the right member is again expanded in a Taylor series and the explicit 
forms of Aq and Ap are inserted, an explicit representation E(q, p, A) is 
obtained. Successive differentiations with respect to p can be used to 
yield the symmetric equation 

3E 3S 1 3 1 2 E 3S 3S 1 3 3 E 3S 3S 3S , 

E(q ’ P ’ X} + ^ + 2 3 Pi 3 Pj 3c k 3c lj + 6 3 Pi 3 Pj 3 Pk 3q i 3q k 


= F(q, p, X) 


3S 3F 3S 1 3 2 F 3S 3S 

3X + 3q i d-p ± + 2 3q t 3q^ 3p i 3p^ 


1 3 3 F 3S 3S 3S 

6 3q t 3 q j 3q k 3p i 3p^ 3p k + ' ‘ ' 

and only the new variables q and p appear. 

Now let n be the mean motion of the body being studied (moon or 
satellite in the present paper) and let m = v/n be a small quantity (i.e., 
m << 1) . Introducing series expansions 
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F 





S 


k 


A q = 2 M 

k=i . 


A P 



v 


with the subscript, k, denoting Ofm^), permits the separation of the equa- 
tions according to order of magnitude, that is, according to powers of m. 


Thus, the explicit transformation equations become 


( A lP) i 


CA 2 pD ± 

C A 3 p)i 


dS i 

3ql 

9S 2 9 2 S ^ 9S j 

+ 9 q i 3 Pj 3 qJ 

3S 3 9 2 s 2 3S! 32 s! 9S 2 

9q7 + 9q i 9pj 9qT + dq^ 9p.. 9qT 

3 2 S! 3 2 Si 9 S 3 
9 Pj 9 P k 3% 

1 3 3 Si 9Si 9Si 

2 9q t 3pj 3p k 9qT 9^ 


and so forth, for Ap, with analogous equations for Aq. 

The equation for the Hamiltonian separates into the "von Zeipel 
equations": 


E o F o 

8E 0 3Sj 9Fq 9^1 

p > 

3E 0 9S 2 3E 1 9S l 1 9 2 E 0 dS 1 dSi 

E 2 + 3p i 9q ± + 9p.j_ 9q i 2 9p i 9pj 9q i 9q_. 

3Si 9F 0 9S 2 9F : 9 S x \ 9 2 F 0 9S], dSi 

= F2 _ v _ + ____ + ____ + _ 9q _ 3q _ ^T- 9^7 
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^ 3E 0 3Sg SE-^ ^^2 ^^2 * 2 *o 3S2 \ 9^Ex 3Sx 3Sx 

E 3 + 9^r 3q7 + 3^ 9 ^T + 9^7 Iq7 + TpTT^T 3q7 3^7 + 2 3p7^P^ 3q^ ^ 

1 3^Eo 3Sx 3Sx 3Sx 
+ 6 3p i dip. 3p k 7q7 3^7 3q^ 


3S 2 

3F 0 3S 3 

3F x 3S 2 

3 F 2 3S| 

3 2 F 0 

3Sx 

3X 

+ 9 q i 9'Pi 

+ 9qi 3 Pi + 

3p t 

3q k 3q^ 

3 Pi 3 Pj 

3 2 F x 

3 S x 3 S ^ 

1 9 3 F q 

3S X 3S X 

3S X 


■i 3c lj 

3 Pi 3p.. 

6 3q i 3q_. 

3q k 3p i 3p.. 

3 Pk 



and so forth. 

It may be remarked that this representation of the transformation 
equations and of the Hamiltonian explicitly in terms of the new variables q 
and p is not the usual practice. Brouwer (ref. 1) uses the mixed set of 
variables q and y. 

In references 1, 2, and 4 the von Zeipel transformation is used succes- 
sively to eliminate only one angle variable, p., at a time. In the applica- 
tions in the present paper, a single transformation will be exhibited that 
eliminates X and all three angle variables simultaneously. Thus, the par- 
tial derivatives disappear from the left members of the von Zeipel equations, 
giving simply the components, E^, explicitly. The equations will be used here 
only in this simplified form. 

These equations can be regarded as a simultaneous set of partial differ- 
ential equations for the components, S^, of the determining function, S. A 
technique that is similar to the classical one of separation of variables will 
be used to construct a recursive set of equations, each of which is a linear, 
first-order partial differential equation with constant coefficients. Speci- 
fically, the angle variables, p, will be the Delaunay variables (7, g, h) , 
respectively, the mean anomaly, argument of perigee, and longitude of the node 
Each term, S^, will be decomposed into the sum of three terms 

S k = Xk + Y k + z k 


where 

X a periodic function of l 

Y a periodic function of g, h, and X, with l absent and the sum, h + X, 

excluded 

Z a periodic function of two variables, g and the sum h + X, with l 
absent 


7 



II 


Similarly, each 
of four terms . 

term, F^, of the Hamiltonian will be decomposed into the sum 

F k = 

with the prefix E denoting a 
angle variables ( Z- , g, h) , the 
above. The calculations will 

XF k + YF k + ZF k + EF k 

function that is independent of X and all the 
other prefixes having the meaning assigned 
then proceed according to the scheme: 

Order 0: 

E o 


Order 1: 

El 

Xl 

Order 2: 

e 2 

X 2 Y i 

Order k > 3: 

E k 

x k Y k-i Z k-2 


This is reminiscent of the schemes of Brouwer and Hori (refs . 1 and 2) 
which yield short-period terms to a higher order than long-period terms . 

Later it will be shown that the components X, Y, Z are obtained as follows: 


The component X k will be obtained by simple quadrature: 


X k = y*(function of type X)d 1 


The component Y k will satisfy an equation of the type 


3Y 9Y 
3A 3h 


f (e) 


6 = i\ + jh + kg, i ? j, i, j, k integers 


Clearly, the solution is 


Y = /f(6)d0 

The component Z k will satisfy an equation of the form 
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az 
n 9h 


f (0) 


9 = ig + j(h + A), i, j integers 


and the solution is 
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z = 


Si + nj 


/ f(0)de 


where £ and n are functions of the action variables (L, G, H) and it will be 
proved that the denominator never vanishes, so that no critical case occurs in 
the lunar theory. In satellite theory, of course, the case of critical 
inclination has to be excluded. 

Lunar theory will be treated first, and then it will be shown that 
satellite theory is simply a degenerate case. 

It is clear from the preceding discussion that the first task is to 
obtain a series representation of the Hamiltonian, in terms of the Delaunay 
variables, and then to exhibit the decomposition by type (X, Y, Z, E) . This 
will be accomplished in the next section. 


LUNAR THEORY 


Development of the Hamiltonian 

The lunar theory to be developed here is that of the motion of the moon 
about the earth, under the assumption that the earth-moon barycenter (B) moves 
about the sun in a Keplerian orbit, the plane of which will be referred to as 
the ecliptic. Earth, moon and sun are regarded as point masses, and planetary 
perturbations are ignored. The Hamiltonian for the classical Delaunay 
variables (L, G, H) and (7, g, ft) is 



F* 



r 2 P 2 (cos ifO 
P J 


+ 


k 2 S E - M 
p4 E + M 


r 3 P 3 (cos . . . 


(see ref. 3, pp. 271 and 291), and 
the law of cosines of spherical 
trigonometry gives (see sketch) : 

cos \p = cos (f + g) cos h 

- sin(f + g)sin h cos I 

where h is the elongation of the 
node, measured from the sun: 


h = ft - (r + 4>) 
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If h is used, rather than ft, as the conjugate variable to H, then the 
Hamiltonian is simply 

F = F' + H 

dt 

since 


dh _ dQ d<J> 3F' 3 / d<J> \ 3F 

dt dt dt = " 3H 3H \dt ' / 3H 

and 

dH _ 3F' 3F' 3F 

dt 3ft 3h 3h 

This choice of variables makes ip independent of the time, so that solar 
effects enter only via p, the distance from the sun to the earth-moon bary- 
center. Many Fourier series expansions and multiplications are thus avoided. 

The usual methods of the theory of elliptic motion (ref. 3, ch. 2, and 
ref. 6, ch. 3) can be used to obtain the solar factors in closed form and in 
series. The only ones needed here are 

'■d * » 

i-f = v 2 (l . 6) 

P 

where v is the sun's mean motion and 


v 2 ’ “ v l' ■ S .T. M • 


0 - eV /2 


6 = — — = yi 6p. cos k A 


dA 


3 = 


d(<j> - A + e sin <j>) 


dA 


00 

£ e k 


cos k A 


k=l 


and 6 k , 6^ are power series starting with e^; for example. 


= 2e - j e 3 . . . 


, 5 2 11 4 

62 = 2 £ " 24 £ • * 
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3l 


= 3c 




To obtain the desired form for the lunar factors, write 
cos = cos (f + g + h) + y sin(f + g)sin h 


where 

1 T *1 H 

Y = 1 - COS I = 1 - 7T 


Then the Legendre polynomials can be written as 

n 

P n (cos <P) = Z P n>k Y k 
k=o 

where p is a trigonometric polynomial in f, with coefficients that are 
n, K 

trigonometric functions of g and h; for example, 

P2 o = T + T c °s(2g + 2h)cos 2f - — sin(2g + 2h)sin 2f 

3 3 3 

P2 , 1 = - 2" sin 2 h + sin h sin(2g + h)cos 2f + -g- sin h cos(2g + h)sin 2f 

3 3 3 

V 2 2 = s ^- n2 ^ ~ 4" h cos 2g cos 2£ + — sin 2 h sin 2g sin 2f 


The Hamiltonian can now be written as a series 


where 


F = 





F 1 = vH 


F, = v6H + v 9 r 2 p 
2 2 *2,0 

F 3 = v 2 3r 2 p 2 ^ 0 + v 2Y r 2 p 2jl 

and so forth, under the assumptions that e = 0(m), y = 0(m). 
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Finally, to separate each term by type (X, Y, Z, E) note that the true 
anomaly, f, enters in the form 

r n cos kf, r n sin kf , 0 <_ k <_ n 

and these can be expressed as trigonometric polynomials in u, the eccentric 
anomaly, the coefficients being functions of a and e (see refs. 2 and 4). 
For example, 

r 2 = a 2 (l + j e 2 ) + A 0 
r 2 cos 2f = ~ a 2 e 2 + B 0 
r 2 sin 2f = C Q 

where 

A Q = a 2 £-2e^- e + cos u^ + y e 2 cos 2uJ 

B 0 = a 2 £-2e0- e + cos u^ + (l - j e 2 ^cos 2uj 

C Q = a 2 /l - e 2 (-2e sin u + sin 2u) 

and each of these is periodic in l, the mean anomaly, with vanishing mean 
value (type X) . 

The desired representation of the Hamiltonian, then, is 

F ° - & Type E 

Fj = vH Type E 

F 2 : xf 2 = j v 2 [A 0 + 3B 0 cos (2g + 2h) - 3C 0 sin(2g + 2h)] 

YF 2 = v6H + v 2 a 2 e 2 cos(2g + 2h) 

O 

zf 2 = 0 

EF 2 = I v 2 a 2 (l + | e 2 ) 
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F3: XF 3 = 3 • XF2 + -j YV2 sin h[-A 0 sin h + B Q sin( 2 g + h) + C Q cos ( 2 g + h)] 

YF3 = yg v 2 3 1 a 2 e 2 [cos ( 2 g + 2 h + X) + cos( 2 g + 2 h - X)] 

+ y yv 2 a 2 ^1 + j e 2 ^cos 2h - e 2 cos(2g + 2h)J + 3 • EF2 

ZF 3 = yy- yv 2 a 2 e 2 cos 2 g 

EF 3 = - Jyv 2 a 2 (l + | e 2 ) 

F 4 : ZF 4 = yg v 2 3 2 a 2 e 2 cos( 2 g + 2h + 2 X) - yg y 2 v 2 a 2 e 2 cos 2g 

EF 4 = | y 2 v 2 a 2 (l+| e 2 ) 

The X and Y components of F4 will not be given, since the analysis will 
not be carried far enough in the present paper to require them. 

It may be remarked that F3 contains the term v 2 3r 2 p 2 Q , which yields, 
among others, the term 

v 2 3a 2 e 2 cos(2g + 2h) 

and here the series for 3 must be used, giving 

v 2 a 2 e 2 3 1 f cos ( 2 g + 2 h + X) + cos( 2 g + 2 h - A)] 

+ li" v 2 a2e2 ^ 2 [ cos C 2 g + 2h + 2X ) + cos(2g + 2h - 2 A ) ] 

+ . . . 

The first line has been included in YF3, and the first term of the second 
line is included in ZF4. There are also higher order contributions that must 
be properly assigned in a complete theory. 

Solution of the Lunar Problem 

The lunar problem can now be solved by means of the generalized 
von Zeipel transformation from the old variables 

x = (L, G, H), y = (l, g, h) 
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to the new variables 


q = (L, G, H), p = (T, g, h) 

where the bar notation is used to emphasize the fact that the new variables 
are the mean values of the old ones, in the usual astronomical sense. That 
is, the difference between each new variable and the corresponding old one is 
a sum of terms each of which is a periodic function of the time with vanishing 
mean value. 


Since the new Hamiltonian contains only the action variables q, the new 
canonical equations are 


d( *i _ 3E_ 
dt 3p^ 


0 


so that the q*s (L, G, H) are constants. Then 


d P; 


dt 


3E 

3q. 


constant 


and the p ! s (Z, g, h) are linear functions of the time. 

Since the von Zeipel equations connecting E, S, and F contain only the 
new variables, it is convenient to omit the bars and write simply (L, G, H) 
and (Z, g, h) during the calculations. Of course, the bars must be restored 
before the explicit transformation equations and the new canonical equations 
can be written. 


It is convenient to divide the discussion into two parts, since the 
equations of order 0, 1, and 2 are degenerate. It is also possible to omit 
many of the algebraic details by recalling that, for example, is the only 

part of Sk that contains Z, and that 


3Fl 8Z k 3Z k 

3H 3h V 3X " 


since = vH and contains h and X only via their sum, h + X. 
Initial stages .- The first von Zeipel equation is simply 


Eo = Po = 


21 / 


and the second is 


3F 0 3X, 1,2 ax. 


3 u 


which separates immediately into Xi = 0, Ex = vH, since Fi is entirely of 
type E. 
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The equation of second order is, in view of the result just obtained 


y 2 3X 

E 2 = ?2 - 7T TT + v 


(»_Ll 

\3h 3A / 


db \3h 3A 

Since ZF 2 = 0, this equation separates into three: 


and integrating gives 


9 x 2 

L 3 


C ^2- = XF 2 


\3X 3h ) 


YF' 


e 2 = ef 2 


Y 1 L 5 

*2 = X v 2 — 
4 P 

[Ai n 

Yi = (+ - A)H 

15 

16 

E 2 = j v 2 a 2 (l 

3 

+ 2 ‘ 


e 2 sin 2u - yy e3 s i n 3u 


where 

A, . /a o d l = /a o C 1 - e cos u)du 
= a 2 £0- e 3 - 2e^sin u + ^ 

Bi = /B 0 dl 

= e 3 - y ejs in u + j e 2 ^sin 2u + e 3 - y e^sin 3uJ 

C x =/c o dZ 

= a 2 /l - e 2 e e + cos u ) ~ y a + e2 ) cos 2u + ^- e cos 3uJ 

Following Hori (ref. 2) the term (1/2) e is added to cos u to 
annihilate the mean value with respect to t. 

The general case, order k >. 3.- When the explicit equation for the new 
Hamiltonian was separated according to the order k, it was implicitly assumed 
that the order of a term is not affected by the process of differentiation. 
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However, the present Hamiltonian, beginning with F 3 contains the factor 
y = 1 - (H/G), so that differentiation may "lose" one power of y, thereby 
reducing the order. For example, the terms 


9F 3 3y 
9y 9G 


and 


3F 3 9y 
9y 9H 


are of second order. When this is taken into account, the general equation, 
of order k can be written in the form 


F k = 


p 2 3X k 


n 


+ V 


rv* 

3 Y k-i\ 

.(*!» 

3F 3 

+ 

9y ^ 

| 9Zk " 2 

V 9h 

ax J 

\3G 

3y 

3G j 

1 9g 


3F, 

9H” 


9F. 

9y 


ll 

3H 


3 Z k- 


3h 


Qk 


where Q k is a sum of terms involving F k and previously computed quantities. 
The currently undetermined components X k , Y] c _ 1 , Z k _ 2 occur onl y in the terms 
exhibited here explicitly. 

Recalling the separation of F 2 and F 3 into components, it can be seen 
that the terms 


/9XF 2 9XF 3 3Z k _ 2 /3XF 2 9XF 3 9 y \9Z k _ 2 

^ 3G + 3y 3G J 9g + ^ 3H + 3y 9H/ 3h 

are entirely of type X, while the terms 


'9YF 2 9YF 3 3 y\ 3Z k _2 /3YF 2 9YF 3 9y \ 3Z k _ 2 

9 G + 9y 3Gy 3g + l 9H + 9y 3H J 9h 


are entirely of type Y. Since ZF 3 contains the factor e 2 as well as 
yv 2 , the terms containing 9 ZF 3 / 3 y and Z k _ 2 can be regarded as of order 
k + 1 , under the assumption that e 2 = 0 (m); that is, these terms are deferred 
to the next stage. Since ZF 2 = (3EF 2 /3H) = 0, the Z component of the 
equation for E k becomes simply 


3EF. 

Icf 


3EF 

~3y 


3_ 3y 
3G 


9Z 


k -2 


3g 


9EF, 


3y 3 ^k-2 


9y 3H 


ZQ k = 0 


Inserting EF 2 and EF 3 gives 


3Z k _. 

3g 


+ n 


9 Z k-; 

9h 


= z Qk 


16 


J 



where 


5 = 


3 a " 1 

4 V2 G 




} z + (1 


Y)(i - § e 2 )] 


3 a 2 /. 3 2 \ 

■ ' 4 U 2 T V 1 * 2 6 ) 


If the typical term of ZQj. has the form f (0 ) with 0 = ig + j(h + X), then 
the corresponding term in Z^_ 2 , is 


f f(6)de 


Si + nj 

The only case that can produce a small denominator is j = 2i, giving 

a * nj = - J V 2 [| e 2 + Y (l + j e 2 )] 

If y vanishes, the orbit plane coincides with the ecliptic, the node loses 
its identity, and g, h, occur only in the sum g + h. Hence, terms of type 
Z must contain the single argument, g + h + A, so that j = i and the small 
divisor does not occur. Hence, the divisor never vanishes. 

Once Zk _2 has been obtained, the equation for Ek separates into three: 


/ 3 Y k-i 

\ ax 


3Y k . 


ah 


87 ~ 

sum 

of 

terms 

of 

type 

")■ 

sum 

of 

terms 

of 

type 

E k = 

sum 

of 

terms 

of 

type 


and these can be integrated by quadrature. 

The equations o f order 3 and 4 .- If the notation of the preceding section 
is used, the known terms of order three are 


/ax 2 

ax 2 \ 

+ ( d 2l 

8F 3 

dy\ 

1 9Yl + 

/3F 2 

3F 3 

a Y \ 

Vah 

ax / 

\3G 

dy 

3G ) 

3g + 

\3H 

3Y 

' 3H / 


First consider the terms involving F 3 and Y lt These can be put in the form 

( Ml 9YA y 9F 3 3Yi 

\dg 3h ) ~ G 


i. 3F 3 / 3Y 1 
G 3y \3g 


ay Sg 
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by noting that 


3y _ _H_ _ 1 - y 
3G q 2 G 


3y _ 1_ 

3H " " G 

The second term, due to the reappearance of the "lost" y, is of fourth order, 
and hence can be deferred to the next stage (this is true at every stage) . 

The first term vanishes, since Yi contains the single argument g + h (this 
does not occur at later stages) . 

Next consider the term 


^TW-= cos (2g ♦ 2h) 


i H 

Y 1 ' G 


H = G(1 - y) 


Again the series expansion for the solar factor 6 must be used, giving 

- V2‘5 1 a 2 e 2 [cos (2g + 2h + X) + cos(2g + 2h - X)] 

- V 2 <S 2 a 2 e 2 [cos (2g + 2h + 2X) + cos(2g + 2h - 2X ) ] 


The terms on the first line are of third order, while the remaining terms 
can be deferred to later stages . The treatment of the other terms in Q 3 is 
straightforward. In particular, the only term of type Z is 


15 


ZQ 3 = ZF 3 = yv 2 a 2 e 2 cos 2g 


and the general theory of the preceding section, with i = 2, j = 0, gives 


Zi = 


yGe 2 sin 2g 


1 * T e2 ‘ \ T ( 1 * I e2 ) 


Expanding the denominator by the binomial theorem, with 


C =iy(l *§e 2 ) 
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gives 

Zi = | yGe 2 sin 2g 

and the remaining terms 

z x (? + + ? 3 + . . .) 

can be assigned to the higher stages Z 2 , Z 3 , . . . . 

The process of inserting Z\ in the terms not used in its calculation is 
straightforward. In particular, the deferred terms involving 9ZF 3 /3y become 


3y dZi _ 75 

By 3G 9g 64 


yv 2 a 2 e 4 (l - y) (1 + cos 4g) 


thus making contributions of types Z and E to the fourth-order equation. 


The rest of the procedure is straightforward, provided that wherever H/G 
occurs it is to be replaced by 1 - y and the terms with the extra y are 
reassigned to the next stage. Then the third-order equation yields 


x 3 


3 L 3 / 

8X 2 + j yv 2 sin h(l 



cos 



sin h 


+ Bi sin(2g + h) + Ci cos(2g + h)] 

- 4 vv 2 [B 2 sin(2g + 2h) + C 2 cos (2g + 2h)] 

* y 

+ 7 'T [7 y ® 2g cos 2g ~ v l &2e2 cos(2g + 2h)J 
v TSAi 3B, 3 C x 1 

LiG 1 + 3 W" COS(2 ^ + 2h ^ - 3 3G 1 sin < 2g + 2h ^J 

Y 2 = f Via 2 6 + ^ e 2 ^ («J> - X + £ sin cj>) 

- ff V! 2 a 4 e 2 |sin(2g + 2h) + j sin(4g + 4h)J 

+ ~ viC^i - P 1 )a 2 e 2 [sin(2g + 2h + X) + -f sin(2g + 2h - X)J 

- g- yv x a z fl - -g- e z + eysm 2h 

+ yvia 2 e 2 (2 - e 2 )sin(2g + 2h) + yvia 2 e 2 ^1 + j e 2 ^sin(2g - 2h) 

+ Yjg yv!a 2 e 2 (2 - e 2 )sin(4g + 2h) 
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where 


c 225 4 2 G 3 2/1 3 2\ 

E 3 = -54- viv 2 a 4 e z -y - 4 YV 2 a 2 (l + j e z ) 

jL 

B2 = 1 d Z and C2 = fc 1 dZ 

Partial results from the fourth stage are 

Z 2 = -| Ge 2 ( 3 2 - 6 2 + 6 X 2 - SiB^sinCZg + 2h + 2X) 

* [- J2 * f ,f2Ge2 ( 1 " sf e2 + ik e ") 

45 _ 9 - 9 ^1 O ^ T . 9 T T 4 •."] - « 


E u = 



256 


[25 

+ 

|_128 

+ 

225 


1024 


125 


3072 

v 2 

2&6 ( 


yGe 


4 - i§§4 Y 2 Ge 2 (16 - 16e 2 + e 4 ) 


yvia 2 e 4 (l - e 2 )J sin 4g 
y 2 e 4 (4 - e 2 )sin 6g 

/ 49 873 2 4347 4 \ 675 2 a 6 e 2 r 2 . 

(- 64 + -sr e - -sir e ) * 1024 v > ”2 T 2- (4 - 5e 5 


9 . 4 , 

512 G 


+ a 4 (16 - 352e 2 + 411e 4 - 75e 6 ) + ^ yv 2 a 2 e 4 


+ y 2 v 2 a 2 (64 - 304e 2 + 400e 4 - 25e 6 ) 


This completes the discussion of the lunar theory. The validity of the 
method has been established, and the calculations have been carried far enough 
to illustrate the manipulative techniques that are required. 


ARTIFICIAL SATELLITE THEORY 


Development of the Hamiltonian 

The primary purpose here is to demonstrate that the generalized 
von Zeipel transformation can be applied directly, in a simple, degenerate 
form, to artificial satellite theory. The usual assumptions about the gravi- 
tational field of a triaxial ellipsoid give the Hamiltonian 
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+ vH + 


F 



CO 

2 >n- 

n=2 


oo n 


EE 

n=2 m=l 


T n,m 


where the zonal harmonics, D n , and the tesseral harmonics, T n#m are 


D n = ‘P J n -5; T P n^ sin 0 ) 
j^n 

T n ,m = P (C njm cos mA + S n>m sin mA)P n m (sin 9) 


Here v is the spin velocity of the earth, R is the equatorial radius of the 
earth, 0 and A are the geographical latitude and longitude (measured east 
from Greenwich), J, C, and S are dimensionless constants, P n is the Legendre 
polynomial, and P R m is the associated Legendre function 


P n m (z) = (l - z 2 ) m/2 


d m P n (Z) 

dZ m 


The notation for the zonal harmonics is Brouwer’s (ref. 1) and the 
notation for the tesseral harmonics is Izsak’s (ref. 7). As in the lunar 
theory the term vH is added to the Hamiltonian because the longitude of the 
node, h, is measured from a rotating reference line (the meridian of 
Greenwich) . 


Orders of magnitude can be defined by introducing n, the mean motion of 
the satellite, the dimensionless parameter m =• v/n << 1, and the usual 
assumptions (Garfinkel, ref. 8) 


Thus , 


J 2 = 0(m 2 ) 

J n = 0(J 2 2 ) , n > 2 

C n ,m = 0(J 2 2 ) 

S nj m = 0(J 2 2 ) 


F = F q + 


F o 

- M 

2L 2 

F 1 

= vH 

f 2 

CSl 

Q 

If 

f 3 

= 0 

F4 

= sum 


and 


Fi + F 2 + F 3 + F 4 


of tesseral harmonics 
remaining zonal harmonics 
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To decompose by types (X, Y, Z, E) 


F k = XF k + YF k + ZF k + EF k 


note that, since the present Hamiltonian does not contain the time, the 
definitions reduce to 


X: 

periodic 

in 

Y: 

periodic 

in 

Z: 

periodic 

in 

E: 

independent 


l 

g and h, independent of t 
g, independent of t and h 
of all three angle variables 



5 2 3 

P 3 (sin 0) = sin 3 0 - sin 0 

= sin 1(1 - 5 cos 2 I)sin(f + g) - sin 3 I sin(3f + 3g) 

P 4 (sin 0) = sin 4 0 - sin 2 0 + |- 

= ^ (3 - 30 cos 2 1+35 cos 4 I) 

- ^ sin 2 1(1-7 cos 2 I) cos (2f + 2g) + sin 4 I cos (4f + 4g) 
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and the associated Legendre functions can be treated similarly, since 


P n m (sin 0) 


cos m 0 


d m P n (sin 0 ) 
d(sin 0 ) m 


Thus, for example, 

P 2 2 (sin 9) = 3 cos 2 0 
P 3 1 (sin 0) = cos 0^y- sin 2 0 - = -|- 

P 3 2 (sin 0 ) = 15 cos 2 0 sin 0 = 15 cos 
P 3 3 (sin 0 ) = 15 cos 3 0 


cos 0[3-5 cos 2 1-5 sin 2 
2 0 sin I sin(f + g) 


I cos ( 2 f + 2 g) ] 


The factor cos m 0 will disappear by cancellation once the longitudinal 
factors have been obtained. To show this, use the trigonometric formulas 
(ref. 9) 


w = A - h 


sin w = 


sin(f + g)cos I 
cos 0 


cos w 


cos (f + g) 
cos 0 


to obtain 


sin A 


ci sin h + si cos h 
cos 0 


. ci cos h - s 1 sin h 

COS A = 

COS 0 


C! = COS (f + g) 


Si = cos I sin(f + g) 


and hence, by induction. 
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c m sin m h + s m cos m h 

sin mA = 

cos m 0 


c m cos m h - s m sin m h 

cos mA = — 

cos m 0 

c m+ l = C m c l - S m s l 
s m+1 = c m s l + s m c l 

so that, for example, 

C 2 = sin 2 I + (1 + cos 2 I)cos(2f + 2g) 

S 2 = cos I sin( 2 f + 2 g) 

C 3 = sin 2 I cos (f + g) + (1 + 3 cos 2 I)cos(3f + 3g) 

3 1 

s 3 = j sin 2 I cos I sin(f + g) + j (3 + cos 2 I)cos I sin(3f + 3g) 

These yield the desired expressions for the tesseral harmonics, typical 
examples being 

T 2 ,2 = |(C 2j 2 cos ^h + S 2> 2 sin 2h) s i ° 2 I + y (1 + cos 2 I)cos(2f + 2g) 

+ (S 2 2 cos 2h - C 2 2 sin 2h)cos I sin( 2 f + 2 g)| 

T 3 , 1 = j h ^ [(C3 } 1 cos h + S 3jl sin h)cos(f + g) 

+ (S 3 y cos h - C 3 j sin h)cos I sin(f + g)] 
x [3 - 5 cos 2 1-5 sin 2 I cos(2f + 2g)] 


Thus, both zonal and tesseral harmonics can be expressed as trigonometric 
polynomials in f + g, multiplied by negative powers of r, and the coeffi- 
cients are functions of I and h. To effect the decomposition by types, 
introduce the mean values with respect to Z (Brouwer, ref. 1): 
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y*27r 

= 2 TT J ° n 

J o 


dl 


-21T 

*n,m 2 tT J T n,m 


dl 


so that the components of type X are simply 


*D n = D n - D n 


XT n,m ~ T n,m T n,m 


Perform the integrations by changing to the true anomaly, f, as the 
variable of integration, using the formulas 


dZ = df , r = 


GL : 


1 + e cos f 


Then the negative powers of r become positive powers of 1 + e cos f, and 
the integrands become trigonometric polynomials in f. 

The zonal harmonics have no components of type Y, and the tesseral 
harmonics have no components of type Z or E. Hence, 

D n = ZD n + ED n 

^n,m ” ^ T n,m 

and the separation of D n is done by inspection. Typical results are: 


F 2 = -p J 2 P 2 C s i n 0) = \ yJ 2 [3 cos 2 1-1+3 sin 2 I cos(2f + 2g)] 

r 3 


4 r 3 


EF 2 = F 2 = ~ J 2 (3 cos 2 I - 1) 

H G°L° 


YF 2 = ZF 2 = 0, XF 2 = F 2 - ef 2 


The third zonal harmonic yields 


n _ 3 T p 5 R 3 e 


ZD 3 = D 3 = ± J 3 


5 t 3 


G b L 


sin 1(5 cos 2 I - l)sin g 
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yd 3 = ed 3 = 0 


xd 3 = d 3 - zd 3 


and the fourth gives 

ZD 4 = J 4 sin 2 1(1-7 cos 2 I)cos 2g 

64 G 7 L 3 

ED 4 = - JL- J 4 (l + % e 2 \ (3 - 30 cos 2 1 + 35 cos 4 I) 

64 G 7 L 3 \ 1 / 

YD 4 = 0, XD 4 = D 4 - (ZD 4 + ED 4 ) 


The first two tesseral harmonics give 

3 . . 4 p 2 

YT 2 2 = sin 2 I (C 2 2 cos 2h + S 2 ^ 2 sin 2h) 

l q 3 l 3 ’ 

XT 2 ,2 = T 2 ,2 - YT 2j2 

and 

5 ^ 

YT 3> 1 = 3- 4 s l 3 e - [ CC3, 1 cos h + S 3j 1 sin h) (1 - 5 cos 2 I)cos g 
+ (S 3 1 cos h - C 3 1 sin h) (11 - 15 cos 2 I)sin g] 

9 > 

It is clear that every harmonic can be decomposed in this fashion. 


Solution of the Artificial Satellite Problem 


Initial stages.- The generalized von Zeipel procedure gives, for the 
first three stages, essentially the first-order results obtained by Brouwer 
(ref. 1). 


E 


o 



Ei 


El 


y 2 

L 3 


= Fi = vH 


and, as in the lunar theory, X x = 0. The second-order equation is 


e 2 


xf 2 + ef 2 


u 2 3X ? 3Y, 

■ + V 3h“ 
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Separating and integrating gives 


Yi = 0 



E 2 


\ ) 


= J J 2 | (3 cos 2 I - 1) (£ - l + e sin f) 

+ 3 sin 2 iT-j sin(2f + 2g) + | e sin(f + 2g) + e sin(3f + 2g) 


The general ca se, order k >. 3 .- Since the quantities X and y of the 
lunar theory do not occur in the satellite problem, the von Zeipel equations 
do not suffer any loss of order due to differentiation. Hence, the general 
equation becomes 


v 2 ax k ^ 

E k = - 7T TT + v 


9 Y k-i 9F 2 9 Z k -2 
3h + 3G 3g 


+ Qv 


where Q k is a sum of known terms. The component Z k _ 2 is obtained from 


3EF. 3Z k 

ZQk - ° 


that is. 


3Z 


3 , y 2 t n 

4 J 2 77 ? (5 cos 1 - l) ~w 


k-2 


= ZQ. 


which can be integrated by quadrature. The general equation then becomes 


Ek ' ' iy n 


2 axi 


3Y k -i 9 *P 2 9 z k- 2 

3h + 3G 3g + 


Qk " 


ZQ k 


which separates into three: 


u 2 9X k 3XF 2 3Z k _„ 

JT TT~ " X< k * Tip ~5g 

3 Y k-i 

V 3h YQ k 
E k = E( ^k 

Notice that the partial differential equations of lunar theory have 
degenerated to ordinary equations here. 



Equations of order 3 and 4.- The third-order equation is completely 
degenerate since the known terms are 


a V 

Q3 = F 3 + V ^ 


3F 2 /3Xj 3YA 3F 2 3X 1 3 2 F D 3Xj 3X £ 

3G - y' 3g~ + 3g ~ ) + 3L~ 3 T~ + ^ 3 7“ = ° 


giving 

E 3 = X 3 = y 2 = Zi = o 

For the fourth-order equation the known terms are 


3F 2 3X 2 3F, 3X 2 


Q 4 F 4 + 9L 3Z + 3G 3g 


1 8fFo /3X 2 V 

2 3L 2 W) 


Since this expression is quite lengthy the components of type X will not be 
exhibited. The zonal harmonics give the components of types Z and E which 
yield Brouwer's results (ref. 1) 


_ 1 y 2 R 2 e 2 sin 2 I r. n , c -n 

22 = 32-^5 cos 2 I - 1 L 2( ' 15 C ° S T) 

+ (1-7 cos 2 I)J sin 2 S ~ j sin 1 cos S 


E 4 = m J 2 


2 H 


6 R 4 


7t 3 


G 7 L 


(5 cos 4 1-18 cos 2 1+5) 


+ 4 £ (3 cos 2 I - l) 2 + 5(7 cos 4 I + 2 cos 2 I - 1) 

+ J 4 ^TT (s - s') (35 cos 4 I - 30 cos 2 I + 3) 

128 G 7 L 3 V l 2 j 

The first two tesseral harmonics yield the typical contribution to Y 3 : 

Y 3 = - f v'g^Li 3 Sin2 I( - C2 » 2 Sin 2h " Sz > 2 C ° S 2h ^ 

+ 1~ y r I [(E COS 2 I - l)cos g(C 3 I sin h - S 3 i cos h) 

° vG b L 3 ’ 

+ (15 cos 2 I - ll)sin g(S 3) i sin h + C 3> i cos h) ] 

Again, the development can be terminated with the remark that the 
validity of the method has been established, provided the critical inclination 
is avoided. 
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CONCLUDING REMARKS 


The generalized von Zeipel transformation presented here is clearly 
capable of producing a highly efficient lunar theory. The major simplifica- 
tion resides in the partial uncoupling of solar and lunar factors. This 
occurs in the initial development of the Hamiltonian and in the final develop- 
ment of the determining function. At both stages many Fourier series expan- 
sions are avoided, the result being a product of two finite expressions. This 
is achieved primarily because the elimination of the short period terms is 
effected by solving an ordinary differential equation, in which solar factors 
play the role of constants. 

If the inclination angle is not restricted to small values, the same 
technique can be applied to eliminate the terms of types X and Y, with terms 
of type Z remaining in the new Hamiltonian. The resulting system might well 
be a fruitful subject for future research. 


Ames Research Center 

National Aeronautics and Space Administration 
Moffett Field, Calif., July 24, 1969 



APPENDIX 


EXPLICIT RECURSIVE ALGORITHMS FOR THE VON ZEIPEL TRANSFORMATION 


The basic, implicit equations of 



the von Zeipel transformation are 


+ 9 S(q, y) 
9/i 

+ 9S (q, y) 
8q. 


where, for present purposes, the dependence on t, via A, is irrelevant. Let 
<f> (x, y) be an arbitrary function and let 

<Kq> p) = <f>(x, y) 

That is, this is an identity in either set of variables under the transforma- 
tion (x, y) -*-* (q, p) . Eliminating p in the left member and x in the 
right member gives 




, r y + y) 
l * L yi sqi 


= 4> 


, + 9S (q, y) 


Cq, y) ~ 

9/i j 


an identity in (q, y) . Expanding each member in Taylor's series gives 


Cq, y) + 


oo 

y± 


^(q. y) 


as as 


as 


k=i 


a yii ay i2 . . . ay ik aq^ aq^ * * ’ a qi 


1 k 


= <f>(q, y) + 


E 

k=i 


as as 


J 3 <Kg, y) 

k! ^ij 9c li 2 • • • 9t }i k 3y ii 9 >% 


as 


3yi 


x k 


with the repeated subscripts ij, i 2> • . satisfying the summation 

convention. Since this is an identity in the variables (q, y) , the variables 
are dummies, and can be replaced by (q, p) or by (x, y) . Introducing the 
series expansion 

oo 

+ = ^n 
n=o 

oo 

^ = Z) ^n 

n=o 
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s = 7 . s. 


and collecting terms of equal order gives, with (q, y) replaced by (q, p) : 


and, for n >_ 1, 


’J’r.Cq, p) = *o(q> p) 


5 k ^.(q, p) 3S a 1 Cq> P) 5S a k (q, P) 


X^k! / j 9p. 3p . ... 9p, 3q.- 

*ii r i 5 ri v ni i 

k=l a 1 2 k i 


= 4> n (q, P) + 


^ k^E 


3 N> a tq, P) SS^Cq, P) 9Sa k (q, p) 

3q. 3q- . . . 3q. 3p n * ' ' ' 3p^ 


where ^ denotes the sum over all combinations of a, a*, a 2 » • * • > a k 


satisfying 


a >. 0 , otj >_ 1 


K 

E 


a + J ^ctj = n 


This can be put in the explicit form 


n " i ^ 

BZE-E 


k=l ai=l a 9 =l 


9pi 3p. ... 3p ± 3q i * Bq, 

J-l i 2 i k 1 1 X 1 


ri r- 


MS-I 


k=l a , =i a 2 = 1 


3q. 3q . . . 3 qi 3p * ‘ * 3 Pi , 
a, =1 12 k l 


where 


ri = n + 1 - k 

r^ + 1 = r j - oij + 1 for j = 1, 2, . . . , k - 1 


a = ri, - a 


k ” k 
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Thus, for example. 


'I' = <f>. 


3 *o 3S l 


= <f> 




3 Pi 9 qi Y 1 3^ Sjk 


3^0 


9S 1 9S! 


9i(j j 9Sj dip 0 9S 2 i 

^2 + 9p^ 9q7 + 9p7 9qT + 2 9p^ 9p ^ 9q7 Bq.. 


= <f>; 


9<f>^ 9Sj 


9 4>o 9S 2 


1 

+ — 


9 2 (() 


9S X 9S X 


9qi 9pi 9q t gpjL 2 dq ± 9q.. 9p t 9pj 


and so forth, a recursive system. Thus, if <f> n Cx , y) are given, then 
ip n (q, P) are obtained recursively. Similarly, if these same equations are 
written with (x, y) in place of (q, p), then d>n C x » y) can be obtained 
recursively when ^ n (q, p) are given. The equations for transforming the 
variables (x, y) and (q, p) , in either direction, now follow immediately as 
special cases, by taking, successively 


<f> 0 ( x , y) 
<?> 0 (x, y) 

(q> p) 

^ 0 (q» p) 

The results are 
where 


X B , 

O 

II 

3 

for 

a 

/ 0 


6 =0 
Y a 

for 

a 

/ 0 

qg> 

Th =0 

Y a 

for 

a 

/ 0 


ip = 0 
a 

for 

a 

/ 0 



q B,o “ q B 


and, for n >_ 2, 

n-l A 

- i!" V' 

q B,n 9 p 3 / v k! jCLj 

k=l ai=l 


q 


6,1 


9Sj 

^3 



a k =1 





9S ai 



9Sa k 

3 ^i k 
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where 


tj = n - k 

tj +1 = tj - ou + 1 for j = 1, 2, . . k - 1 

a = t k - a k + 1 

Thus, for example, 

_ 9S 2 9c lB, 1 3S! 

q 6,2 " 3 p 6 " 3 p t 3 q t 

1 92q e,i 3S 1 9S i 

q 6,3 3p g " 3p i 3Pi 3q t " 2 3 Pi Spj 3q i 3q.. 

and so forth, where = S^Cq, p) . 

The transformation equation for y is 


3 


oo 

= p g,n 
n=o * 


where 


p e,o p g 


0,1 


3 Sx 

^7 


and, for n >. 2, 


0,n 


k=i 

Thus, for example. 


n-i x i k 

BE E 

k=l ai= 1 a i=l 


3 k p 


0,a 


3 S n 


3 Pi 

a k =1 1 


■ 3p iv 9q i, 


3S 


a k 


3 q 


1 k 



3S 2 

9p 3,l 

3S 1 



p 3,2 

9q 3 " 

3p i 

3q. 




3 s 3 

9p 3, 1 

3S 2 

9p e,2 3Si 

i 92p 3,i asj 3S, 

p 3,3 = 

9q 3 

3 Pi 

3q. 

3pi 3q t 

2 3 Pi Sp^ dq i 3qj 

th, again 

with 

S k - 

s k Cq, 

P) • 



The reverse transformation is 


q 3 


- E 

r\—r\ 


X 


3,n 


where 


X 3,o X 3 


3S] 

W f 


and, for n >_ 2, 


3,n 


n-1 r l r k 

£ 


a k x 


3, a 


3S 




3S, 


a k 


3x. 


k=l ai=l 
Thus, for example. 


a, =1 

k 


il 


• 9 *i k 8 n 1 ' ' ' s n, 


P.2 


3S 2 3x 6> , 9Sj 

3 y e - iy L 


8x g, i »S 2 3X 6.2 PSt 1 82x B , i 8S t 3S, 

X 3,3 “ 3y g dx ± By i dx ± dy ± 2 dx ± dx^ 3 yi 3y.. 

and so forth, where now S k = S k (x, y) . 

Finally, the transformation equation for p is 


where 


■S'- 


y 3,o y 3 
3S] 

y 3,i = 93T 


and, for n >_ 2, 
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as. 

r 3,n ~ 3x~ 


n 


n-i "1 k 

2>E---E 

k=i ai=i a k =l 


a k y 


3, a 


as, 


“i 


ax, 


^ • • • 3 x i k 3y,^ 


as, 


a k 


3y ir 


Thus, for example. 


as, a y 


3,1 9S 1 


} 3,2 ax g axi ay, 

_ as 3 9y B ,l 3S 2 3y 3j2 as, ! 3 %,1 3Sj 3Si 

y 3,3 ax 3 ' ax, ay, “ ax t ay, ' 2 ax, ax.. ay, ay^ 


and so forth, again with = S^fx, y) 


It may be noted that the authors of recent papers on the Lie transforma- 
tion (refs. 10 and 11) have remarked that the von Zeipel transformation is 
unsatisfactory because of its implicit nature. Such remarks are no longer 
valid in view of the present explicit algorithms for effecting the transforma- 
tion of coordinates and of arbitrary functions in either direction. 
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